The x- dependence of Parton Distributions Compared with Neutrino Data 
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We use the variational principle of Quantum HadronDynamics, an alternative formulation of 
Quantum ChromoDynamics, to determine the wavefunction of valence quarks in a baryon at a low 
value of Q^. This can be used to predict the structure function xFz{x,Q^) at higher values of 
using the evolution equations of perturbative QCD. This prediction is compared to the measurements 
- - - of neutrino scattering cross-section by the CDHS and CCFR experiments. The agreement is quite 

G^ 

^\ ' good, confirming the validity of QHD as a way of studying hadronic structure 
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One of us has proposed [g| that it is possible to describe strong interactions directly in terms 
of hadrons rather than in terms of quarks and gluons; this new approach was called Quantum 
' HadronDynamics (QHD). It is equivalent to Quantum ChromoDynamics (QCD, the universally 

(N ■ 

QQ ' accepted theory of strong interactions), except that the semi-classical approximation of QHD is a 

o ■ 

0^ ' applied directly to QCD give results that are too inaccurate except in the limit of short distances. 



QHD has been constructed so far only in two space-time dimensions. However, that is sufficient to 
obtain the structure functions of Deep Inelastic Scattering: in the limit of zero transverse momentum 
^ for the constituents, we can dimensionally reduce the theory to two space-time dimensions. More 

precisely, the dependence of the structure functions on the Bjorken variable is determined by ignoring 
the transverse dimensions Thus we can predict the non-perturbative xb dependence of structure 
functions, something that is completely inaccessible to the standard perturbative formulation of 
QCD. The leading order effects of transverse momenta are described by the DGLAP |^ equations 
that determine the dependence of the structure functions: this is the fundamental insight of the 
paper of Altarelli and Parisi. 

In this paper we will focus on understanding the distribution of valence quarks in the nucleon. 
This is measured directly in neutrino (and anti-neutrino scattering) against the nucleon: the isospin 
averaged valence parton distribution function is just the structure function F3{xb) measured in this 
process Q. To leading order in jog^g^) ; it is not necessary to know the gluon structure functions of 
the hadron in the DGLAP evolution of this quantity. The anti-quark content of the baryon is zero at 
the initial value — Qq within our approximations. This is consistent with the phenomenological 
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model of Gliick and Reya [g], that the anti-quark distribution of the proton is dynamically generated 
by evolution from an initial value of zero at low . We will, in a later paper [Q, calculate the 
anti-quark distribution functions and show that this is indeed justified. However, we expect the 
"primordial" gluon distribution to be non-zero. 

In a previous paper we derived a variational principle that determines this valence quark 
distribution function in the semi-classical approximation of QHD. It was also shown there how to 
derive this variational principle from QCD. In a separate paper |^ we derived the same principle 
from an interacting valence parton model of the baryon. In this paper we will obtain approximate 
solutions of the variational principle and compare with the experimental data of the CCFR and 
CDHS 10 collaborations. The agreement of the predictions of QHD with data is quite spectacular: 
we now have a theory of the xb dependence of deep inelastic structure functions. 

Let V'(p) be the wavefunction of a valence parton in a proton, thought of as a function of the null 
component of momentum (p = — p^) j^. For kinematical reasons, P > p > 0, where P is the 
total momentum of the proton. Naively, we would have the sum rules. 



|^(p)|^^ = l, p|^(p)|-_^ = P. 

However, we have to modify the momentum sum rule since it is known that the valence partons 
carry only about half the momentum of the proton: the rest is carried mostly by the gluons, with 
a small contribution from anti-quarks. Thus we impose -^c /p^p|V'(p)P|f ~ fP, where / is the 
fraction of the momentum carried by the valence partons. We will address the problem of gluon 
distribution functions in a later paper, where we will attempt to calculate /; here we will treat it as 
a parameter. 

The wavefunction is determined by minimizing the energy 

subject to the above conditions. (Here i'{x) — ''Pip)^^''^ ^ is the wave-function in position 
space.) This variational principle was derived from QCD through QHD in Ref. as well as from 
an interacting parton model in Ref. It describes partons (within the mean field approximation) 
which are interacting with each other through a linear Coulomb potential, which binds them into 
a baryon. The linear potential comes from eliminating the gluon fields: it is the one-dimensional 
Fourier transform of the gluon propagator. 

A Lorentz invariant formulation (which is more convenient for our purposes) would be to minimize 
the mass^ of the baryon rather than its energy: 



^!^(p)r^ + Y J m^^rmvr^^^dxdy 



We are ignoring the flavor and spin quantum numbers of the partons so what we get will be the 
spin and flavor averaged wavefunction. Also /i^ — — ^ is the quark mass^ after a finite 

renormalization and m is the current quark mass. The dimensional parameter g ^ A.qcd 

determines the strength of the interaction. We now define the number density of valence partons 

i2\ „ /,n2\ „ /'/n2\ 



V{xb) = iVc 



2 



We have normalized this so that the Gross-Lewellyn-Smith sum rule (including the perturbative 
corrections up to order a^{Qo)) is satisfied. The coefficients Ci,C2 and C3 are given in Ref. 
We will finally set the number of colors Nc to 3 and the number of flavours to 2 at the initial low 
value of Qo- 

Since the total momentum scales like P ~ A'^c, in the limit as Nc 00, the parton momentum 
has the range < p < 00. In the limit when Nc 00 and m = 0, we have found an exact minimum 
of the variational principle: ip^p) = Ce"s . The condition for minimizing is an integral equation 
and we can verify that this function is a solution by explicit computation. (The calculation involves 
infrared singular integrals, which are defined through an appropriate principal value prescription 
as in [^.) This suggests that even when Nc is finite, a reasonable variational ansatz would be 
'^(p) = C [1 - f (RecaU that e"^ = lim„_oo[l - f]".) The constant C is determined by the 



normalization condition. By using the momentum sum rule, we get 6=^ — 1 + ol[^ — 1]. "a" is 



determined by the variational principle. It is zero in the limit of chiral symmetry and rises like a 
2 

power of In the limit A^'c ^ 00, a is determined by the transcendental equation 
.™2 /-i 



9^ 



= 1 + 



dy 



(l+j/r + (l-y)"-2 



dy 

y2 



which we derived by a Frobenius-type analysis of the integral equation for the minimization of 
in an earlier paper [^. Thus in the limit m = 0, a = Q, and we have the valence parton 
distribution function V{xb) = C[l — a;6]~r where C is fixed by the GLS sum rule. This variational 
approximation agrees well with our numerical solution of the same problem in Ref. Q but is much 
simpler to use. The limit of chiral symmetry is a good approximation in the case of the nucleon 
since the current quark masses of the up and down quarks (5-8 MeV) are small compared to the 

energy scale of the strong interactions {Aqcd ~ 200 MeV). The valence parton distribution depends 
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on g only through the ratio ^ and in the limit m = 0, is independent of g. 

The structure function Fs{xb, Q^) of neutrino scattering on an isoscalar target has been accurately 
measured by the CCFR and CDHS experimental groups at several values of xb and Q^. To 
compare with that data we need to evolve our computed distribution function to the appropriate 
values of by the DGLAP equation: 

dV{xB,t) as{t) dy xs 



dt 2tv y ' y 

where t = log(Q^/(5o)) Qs(^) is evaluated using the two-loop (3 function |^ and Pqq is the evolution 
kernel to leading order given in Ref. [^]. We solve the evolution equation numerically, assuming an 
initial value of Qg — 0.4 GeV^. This low value of Qq is justified since we can show Q that the 
anti-quark distribution function of the nucleon is quite small. The range of Q over which we are 
evolving is small {Q ~ 0.6 to 5 GeV) and its effect is small, thus justifying the use of the leading 
order DGLAP equation. We also set Nc = 3, Aqcd = 200 MeV and the current quark mass m — 0. 
As for the fraction of baryon momentum carried by valence partons /, we choose it to be 0.5, which 
is in agreement with the phenomenological fits of parton distribution functions ( See e.g., [^.) In 
a later paper we will derive the gluon structure function as well, and at that time we will have a 
theoretical prediction for this parameter. The plots show that we have quite good agreement with 
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data. Thus we have estabhshed that Quantum HadronDynamics is a successful way of deriving 
hadronic structure functions from QCD. 
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Figure 1: Comparison of QHD prediction of xF^ (solid curve) with measurements by CCFR(*) 
and CDHS(^). We have chosen the parameters Ql = 0.4 GeV^ and / = 0.5. (a) CCFR at 
= 7.9 GeV•^ CDHS at 7.13 < < 8.46 GeV^ and QHD prediction at = 7.9 GeV^ 
(b) CCFR at Q2 = 12.6 GeV^, CDHS at 12.05 < < 14.3 GeV^ q^LT) prediction at 

Q2 = 13 GeV'^. 
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Figure 2: Comparison of QHD prediction of xFz (solid curve) with CCFR(») measurements at 
(a) Q2 = 20 GeV'^ and (b) = 31 5 GeV"^ . 
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